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У статті запропоновано метод бінаризації напівтонових зображень з використанням оптимального 
порогу на основі застосування кореляційної функції. За ефективністю метод не поступається відомим 
методам глобальної бінаризації, зокрема, методу Отсу. 
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Кеу УУ/огаз: ітаєе Біпагігайоп, орійта! Фгезроїй, соггеіайоп ГипсПоп. 


В статье рассмотрен метод бинаризации полутоновьх зображений с использованиєм оптимального 
порога на основе применения корреляционной функции. По зффективности метод не уступаєт 
известньм методам глобальной бинаризации, в частности, методу Отсу. 

Ключевь(е слова: бинаризация изображений, оптимальньій порог, корреляционная функция. 


Актуальність теми та постановка задачі 


Велика кількість різноманітних металевих конструкцій та виробів у машино- 
будуванні, літакобудуванні, хімічній та інших галузях вже вичерпали свій ресурс або 
підходять до його межі. Створюються нові матеріали, постійно зростають вимоги до 
їх якості. І в першому, ів другому випадку є необхідність контролювати стан матеріалу, 
тобто його механічні властивості. У першому випадку контроль дає відповідь на запи- 
тання - чи безпечною є подальша експлуатація виробу, у другому випадку - чи відпо- 
відає продукція поставленим вимогам щодо якості. 

Одним із способів оцінки стану та якості матеріалу у цих двох випадках є ви- 
значення вмісту структурних і фазових складових, їх геометричних розмірів, кількості, 
взаємного розташування. Співвідношення, які пов'язують параметри мікроструктури 
сплаву з його фізичними та механічними властивостями в кількісній формі, формують 
простий 1 зрозумілий метод контролю цього матеріалу. Проте проводити такі до- 
слідження на рівні мікроструктури візуально та вручну є дуже складно, а часом 1 не- 
можливо, тому це призводить до необхідності розробки та використання комп'ютери- 
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зованих систем обробки й аналізу металографічних зображень для отримання числової 
оцінки взаємозв'язку «хімічний склад -структура - властивості». До найбільш відомих 
програмних продуктів, які розв'язують питання даного напрямку, належать системи, які 
розроблені фірмами 51АМ5 (1, РАСЕ Тесіпоїоєіез (21, ВідеоТест |3| та ін. У більшості 
з цих систем під час обробки та аналізу зображень є присутнім етап, під час якого 
вхідне напівтонове чи кольорове металографічне зображення перетворюють у бінарне і 
на основі цього визначають геометрію шуканих об'єктів, фрактальну розмірність тощо. 
Саме цей етап, в більшості випадків, є найбільш відповідальним у дослідженні мікро- 
структури матеріалів, оскільки від коректності його проведення залежить точність отри- 
маних результатів в цілому. Дана стаття присвячена питанню визначення оптимального 
порогу при перетворенні напівтонового зображення у бінарне. 


Огляд відомих підходів 
При бінаризації зображення інтенсивність кожного піксела 1, і) порівнюється 


з пороговим значенням інтенсивності /.,., яке може бути постійним для всього зобра- 


ження чи адаптивним. Якщо значення інтенсивності піксела більше за значення порогу, 
тоді на бінарному зображенні відповідний піксел буде представлений «білим», в іншому 
випадку - «чорним». При такому перетворенні важливим є правильний вибір порого- 
вого рівня, оскільки він суттєво впливає на формування бінарного зображення, а отже 
і на отримані результати його наступного аналізу. На рис. 1 наведено результати (б, 
в та г) обробки вхідного зображення (а) з використанням різних порогових значень. 
Отримані внаслідок обробки зображення рис. 16, в, г є різними, такими ж різними 
будуть 1 результати їх аналізу. Для прикладу, з такими проблемами доводиться зустріча- 
тися при фрактальному аналізі металографічних та фрактографічних зображень, коли 
для обчислення фрактальних розмірностей необхідно використовувати бінарні зобра- 
ження (|4|; подібні проблеми виникають і при вимірюванні розмірів зерен (|5|, де їх 
межі представляються у бінарному вигляді тощо. Це ще раз доводить актуальність 
задачі вибору оптимального порогового рівня при бінаризації зображень. 

Загалом усі методи бінаризації за способом побудови порогової поверхні І, поді- 
ляються на дві групи: методи локальної та глобальної бінаризації. Порогова поверхня 
І, представляє собою матрицю, розміри якої співпадають з розмірами оброблюваного 
зображення. Кожна комірка матриці І, задає поріг бінаризації для відповідної комірки 
на оброблюваному зображенні. В методах глобальної бінаризації порогова поверхня 
утворена постійними значеннями порогової інтенсивності, у методах локальної бінари- 
зації значення порогової поверхні змінюється від піксела до піксела і залежить від ха- 
рактеристик локального околу оброблюваного зображення. 

Розглянемо деякі найбільш вживані методи бінаризації напівтонових зображень. 
У даній роботі зупинимося на глобальних методах, які володіють більшою швидкодією 
у порівнянні з локальними, що є важливим при обробці металографічних та фрактогра- 
фічних зображень, які часто мають доволі великі розміри. 


Рисунок 1 - Обробка вхідного зображення (а) глобальним пороговим методом 


при значенні порогу І, - 0,3 (6), І, - 0,5 (в), І, «0,7 (г) 
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Найпростіший з методів порогової обробки полягає в поділі гістограми зобра- 
ження на дві частини за допомогою єдиного глобального порогу І. : 


бен ОГО А о р 
убтеног а) 
І, І 

де Щі, )) - рівень інтенсивності вхідного зображення І, в точці з координатами 
(і, 7); І, - порогове значення інтенсивності; /,, (і, /) - рівень інтенсивності оброб- 
леного зображення /, в точці з координатами (ії, )). 

Сегментація зображення згідно з цим методом здійснюється шляхом поелемент- 
ного сканування зображення, при цьому кожний піксел відноситься до об'єкта або 
фону, залежно від того, чи перевищує інтенсивність даного піксела значення порогу 
І, чині. Визначення величини порогу за допомогою гістограми інтенсивностей є прос- 


тим та ефективним методом, якщо гістограма зображення носить чітко виражений 
бімодальний характер. Можливі також модифікації розглянутого методу - дворівневі 
та багаторівневі порогові перетворення. Але, на жаль, в більшості випадків металогра- 
фічні та фактографічні зображення мають не бімодальну, а більш складну гістограму, 
тому використання описаного методу не призводить до бажаного результату обробки. 

Іншим широковживаним методом є метод порогової обробки на основі гра- 
дієнта інтенсивності |б6). Алгоритм обчислення порогового значення складається з 
наступних кроків. Спочатку обчислюємо модуль градієнта інтенсивності для кожно- 
го піксела зображення 


Сі, б є тах (6, ЛЬ, ДУ (2) 
де 1, (і, )) - Цін 1, ))- Ці - 1, )), 2 (і, )) є Ці, ук 1)- ЦЕ, ) - 1). 


Далі обчислюємо порогове значення у відповідності з виразом 


М -1М-1 М -1М-1 


г, -| У. У І, СС, 5) 7/1 У У Єтуп) |, 3) 


із у-ї і-ї у-І 


де М,М - розміри вхідного зображення. 
Результат роботи методу представлений на рис. 2. 


Рисунок 2 - Результат застосування методу порогової обробки на основі градієнта 
інтенсивності: а - вхідне зображення, б - оброблене розглянутим методом 


Одним із найбільш уживаних методів глобальної бінаризації є метод Отсу |7. 
Даний метод для пошуку порогу бінаризації використовує гістограму розподілу значень 
інтенсивностей зображення. Порогове значення згідно з цим методом обчислюється 
таким чином, щоб мінімізувати середню помилку бінаризації, тобто середню помилку 
від прийняття рішення про належність пікселів зображення до об'єкта або до фону. 
Таким чином, дане порогове значення є оптимальним з точки зору мінімуму помилки 
бінаризації (рис. 3). 
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Рисунок 3 - Обробка вхідного зображення (а) методом Отсу (6) 


Метод бінаризації з використанням оптимального порогу 


На підставі проведеного аналізу існуючих методів порогової обробки можна 
зробити висновок, що бінаризація є ефективною тоді, коли оброблене зображення 
відображатиме всі деталі вхідного напівтонового зображення. Іншими словами, вхідне 
напівтонове та результуюче бінарне зображення мають бути максимально подібними. 
Саме ця ідея є основою розробленого методу. Для цього вхідне напівтонове зобра- 
ження перетворюють у бінарне згідно з виразом (1), використовуючи всі можливі 
значення порогу І, із заданого діапазону. Далі для кількісної оцінки подібності вхід- 
ного напівтонового та бінарного зображень використаємо коефіцієнт кореляції К: 


У У, 5-2 (ув То) 


желе і (4) 


У У (ше р-1 М 55 /ць 0-1) 


із узі ізї у-і 


де Ці, )), І, (і,.)) - інтенсивності елементів вхідного та бінарного зображень 


у точці з координатами (і, 3) відповідно; І, І,, - середньоарифметичні значення інтен- 
сивностей вхідного та бінарного зображень, які обчислюються за виразами відповідно 


ре МУ, Їх Із оо У бі) (5) 


В результаті проведених обчислень отримаємо залежність між коефіцієнтами 
кореляції (4) та відповідними пороговими рівнями І. (рис. 4). Зауважимо, що дана 
залежність отримана для зображення, яке представлене на рис. Іа. Експериментально 
встановлено, що для інших зображень вигляд графіку отриманої залежності буде 
подібний за формою до отриманого. 

Таким чином, із сказано вище та рис. 4 слідує, що оптимальним у методі біна- 
ризації буде таке порогове значення І, , яке забезпечує максимум коефіцієнта кореляції 
Кк у виразі (4), тобто максимальну подібність вхідного напівтонового та бінарного зобра- 
жень. Результати комп'ютерної реалізації розробленого методу представлені на рис. 5. 
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Рисунок 4 - Залежність коефіцієнта кореляції К від порогового рівня І 


Експериментально встановлено, що за ефективністю обробки запропонований 
алгоритм не поступається відомим методам глобальної бінаризації, зокрема, методу 
Отсу |71. 


Рисунок 5 - Метод бінаризації зображення з використанням кореляційної функції: 
а - вхідне напівтонове зображення; б - результат застосування запропонованого методу 


Висновки 


У даній роботі запропоновано метод бінаризації напівтонових зображень з вико- 
ристанням оптимального порогу, який обчислюється із застосуванням кореляційної 
функції. За ефективністю метод не поступається відомим методам глобальної бінари- 
зації, зокрема, методу Отсу. При цьому слід зазначити, що якість роботи усіх методів 
глобальної бінаризації, у тому числі і запропонованого, суттєво залежить від вхідного 
зображення, зокрема, від рівномірності фону. Тому в деяких випадках, перед викори- 
станням цих методів необхідно спершу провести вирівнювання фону вхідного зобра- 
ження. Одна з переваг запропонованого методу полягає в тому, що тепер на етапі 
бінарної обробки усунуто суб'єктивний фактор, що дозволяє порівнювати результати 
аналізу металографічних зображень, отриманих, зокрема, за різних умов освітлення. 
Запропонований у роботі підхід буде корисним при аналізі металографічних та факто- 
графічних зображень, при обчисленні їх фрактальних розмірностей, виділенні країв, 
геометричному аналізі меж зерен тощо. 
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РОМ ОНА 


КЕ5СМЕ 


І.М. УПигауе! 
Віпагігатоп Метоа ої МаішаПо»єтарпіс Ітаєе5з уд Оріїта! 
Тиге5Поїа 


І позі са5е8, опе ої Ше піапдааогу зкеря іп Ше апа!узбі5 ої дїота! ітаєєз 15 Шеїг 
ірге5роїд ргосезвіпо. ТПі5 таКез їс роз5ібіе (о 5еЇесі їегі5 ОЇ іпіеге5і Їог Гигібег 
ргосезвіпе апа апаїузія. 

Іп ія са5е, Фе обіесі ої іпіегезі аге 5ігисіига! согпропепіз ої плегаЙоєгарріс ітаєез. 
Опе мау (о а85е55 Ше 58кагиз апа дцаПсу ої паїегіа! 15 дегегпліпайоп ої 5їгискига! апа ррба5е 
сотропепіз5 апа Шеїг єеотеїгіса 5175, питрег, геіайуе розійоп. Маше раї Ббіпа Ше 
рагатегегя ої Фе писгобігисіиге ої Ше аПоу ул їі5 рПузіса! апа песпапісаї ргорегиез іп 
даапікайуе Їогт, Їогті а 8ігаріе апа іпішійуе сопігої! тефодй ої Фіз пагегіа!. Егот Ше 
абоуе ГоПом/ 8 Шас Фе ргоседиге їог сопуегіпє апу раїопе шпаєе іп Біпагу 15 уегу 
тезропзібіе. Сотгесіпе55 ої 1і5 бієпійсапі ітрасі оп Ше ассигасу ої Гигірег апаїЇузі8. ТПаї 
Шрі8 1855це 15 деуоіед (0 Шіз8 могК. 

Тре рарег ргезепі5 Ше ш5е ої орйта! Шге5роїд біпагігайоп. ТПе сгі(егіоп ОЇ 
оритайу 15 Фе пахітит сопеїайоп їипсйоп ої Ше іприє раїйопе апа біпагу ітаєе5 
ргосез5ед. Сопдайсіедй сотриїег 5іпаціайопя ої Ше ргороз5ед арргоасіп 5ромуеад ейісіепі 
ргосез8іпе, а|єогійт доез пог уїеід а умеї-Кпом/п 8Їобаї Біпагігайоп плеїродз, іпсіидіпеє 
Оїви тешоа. 
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